Classification of Abelian and Non-Abelian Multilayer 
Fractional Quantum Hall States Through the Pattern of Zeros 
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A large class of fractional quantum Hall (FQH) states can be classified according to their pattern 
of zeros, which describes the way ideal ground state wave functions go to zero as various clusters 
of electrons are brought together. In this paper we generalize this approach to classify multilayer 
FQH states. Such a classification leads to the construction of a class of non- Abelian multilayer FQH 
states that are closely related to gk parafermion conformal field theories, where gk is an affine simple 
Lie algebra. We discuss the possibility of some of the simplest of these non- Abelian states occuring 
in experiments on bilayer FQH systems a.t v = 2/3, 4/5, 4/7, etc. 



I. INTRODUCTION 

One of the most important questions in condensed 
matter physics relates to how we should characterize and 
classify the many different quantum phases of matter. A 
large part of the story is the theory of symmetry break- 
ing, which tells us that we should classify various phases 
based on the symmetries of the ground state wave func- 
tion. Yet with the discovery of the fractional quantum 
Hall (FQH) states"'^'^ came also the understanding that 
there are many distinct and fascinating quantum phases 
of matter, called topological phases, whose characteriza- 
tion has nothing at all to do with symmetry. How should 
we characterize and systematically classify the different 
possible topological phases that may occur in a FQH sys- 
tem? 

Let us run through a few obvious possibilities. We 
know that the FQH states contain topology-dependent 
and topologically stable degenerate ground states, which 
allow us to introduce the concept of topological order 
in FQH states. •^'^ Such topology-dependent degenerate 
ground states suggest that the low energy theories de- 
scribing the FQH states are topological gauge field theo- 
ries in 2-1-1 dimensions, called Chern-Simons theories. 
So one possibility is that we may try to classify the differ- 
ent FQH phases by classifying all of the different possible 
Chern-Simons theories. But this is not a satisfactory ap- 
proach for non- Abelian FQH states because we do not 
have a good way of knowing which Chern-Simons theo- 
ries can possibly correspond to a physical system made 
of electrons and which cannot. 

Another possibility comes through the connection be- 
tween FQH wavefunctions and conformal field theory 
(CFT). It was discovered around 1990 that correlation 
functions in certain two-dimensional conformal field the- 
ories may serve as good model wavefunctions for FQH 
states.^ Thus perhaps we may classify FQH states by 
classifying all of the different conformal field theories that 
may be used to construct a "valid" FQH wavefunction. 
While the connection between FQH wave functions and 
CFT correlators has been extremely fruitful in both con- 
structing new model wavefunctions and understanding 
their topological properties, the CFT approach is incom- 
plete in that there does not exist a complete classification 



of all possible conformal field theories that may be used to 
construct FQH wavefunctions. This precludes the possi- 
bility of systematically classifying FQH states. The CFT 
approach also obscures the essential physics and the es- 
sential properties of the conformal field theory that allow 
one to obtain amenable FQH wavefunctions. 

In an attempt to obtain a systematic classification of 
FQH states, it was shown recently that a wide class 
of single-component FQH states and their topological 
excitations can be classified by their pattern of zeros, 
which describe the way ideal FQH wavefunctions go 
to zero when various clusters of particles are brought 
together. ^""^^ This analysis led to the discovery of some 
new non- Abelian FQH states whose corresponding CFT 
has not yet been identified. It also helped elucidate the 
role of CFT in constructing FQH wavefunctions. The 
CFT encodes the way the wavefunction goes to zero as 
various clusters of electrons are brought together. The 
order of these zeros must satisfy certain conditions and 
the solutions to these conditions correspond to particu- 
lar CFTs. Thus in classifying and characterizing FQH 
states, one can bypass the CFT altogether and proceed 
directly to classifying the different allowed pattern of ze- 
ros and subsequently obtaining the topological properties 
of the quasiparticles from the pattern of zeros. This con- 
struction can then even be thought of as a classification 
of the allowed CFTs that can be used to construct FQH 
states. Furthermore, these considerations give way to 
a natural notion of which pattern of zeros solutions are 
simpler than other ones. In this sense, then, one can 
see that the Moore- Read Pfaffian quantum Hall state^ is 
the "simplest" non- Abelian generalization of the Laugh- 
lin state. 

In this paper, we generalize the pattern-of-zeros clas- 
sification to multilayer FQH wave functions. Such a sys- 
tematic classification leads us to the construction of non- 
Abelian multilayer FQH states and also helps identify the 
simplest non-Abelian generalizations of the Halperin-'^^ 
bilayer FQH states. Therefore, in the same way that the 
Pfaffian FQH state is the simplest non-Abelian single- 
layer state and so one of the first non-Abelian states ex- 
pected to be realized experimentally, our construction al- 
lows us to identify the simplest non-Abelian bilayer states 
and therefore some of the ones that deserve further con- 
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sideration in future experimental and numerical work. 

We would like to point out that the "ID charge- 
density- wave" characterization of single-component FQH 
states"'^'''"'^^ is closely related to the pattern-of-zeros ap- 
proach. Our pattern-of-zeros approach for multilayer 
FQH states can also be viewed as a way to gciiicralize 
the "ID charge-density-wave" characterization to multi- 
component cases. 

This paper is organized as follows. We begin by de- 
scribing the ideal FQH wave functions that we can char- 
acterize by the pattern of zeros and their relation to sym- 
metric holomorphic polynomials. In Section IV, we de- 
scribe the different ways to characterize the pattern of 
zeros. In Section V we find the conditions that the pat- 
tern of zeros must satisfy in order to describe valid FQH 
wave functions. In Section VI, we sketch how one may be- 
gin to construct ideal Hamiltonians whose ground states 
will be FQH wave functions with a given pattern of zeros. 
After a brief summary of the pattern of zeros data and 
conditions in Section VII, wo explain in Section VIII the 
relation between the pattern-of-zeros approach and the 
CFT approach to FQH wave functions. In Section IX, 
we describe some example solutions of this systematic 
classification of multilayer FQH wave functions, which 
yields many non-Abelian multilayer states. In Section 
X we discuss some of the simplest of these non-Abelian 
bilayer states that may be relevant for experiments on 
two-component quantum Hall systems and that warrant 
further numerical study. 



II. FRACTIONAL QUANTUM HALL STATES 
AND SYMMETRIC POLYNOMIALS 

The ground state wave function of a two-dimensional 
system of electrons in the lowest Landau level can be 
written in the form 

* = $(zi,--- ,ZArJe-3J:r=il-^l', (1) 

where Zi = Xi+ iyi, {xi, yi) are the coordinates of the i^^ 
electron, and • • • , Zpj) is a holomorphic function of 

Zi. Since the electrons obey Fermi statistics, <& is anti- 
symmetric under interchange of any two coordinates Zi 
and Zj when all of the electrons are identical. In many 
physical situations, the electrons may be distinguished by 
various quantum numbers, such as a spin index (when 
the Zeeman energy is not too high), a layer index (in 
a multilayer two-dimensional electron system), a valley 
index (such as in graphene and SiGe heterostructures), 
etc . In such cases, the ground state wave function in the 
lowest Landau level may instead be written in the form 

^r = $({4})e-i^*,/l^'l', (2) 

where 7 = 1, • • • ,Nf is a. flavor index and Nf is the num- 
ber of different flavors. $ is then necessarily antisym- 
metric only under interchange of zf and zj for any i and 



j. Given any antisymmetric polynomial 4'anti-sym({2/}), 
we can uniquely construct a symmetric polynomial: 

ff. — ^anti-sym({^/}) 

'i'symm wiU also be a polynomial because $anti-sym must 
vanish when any two identical particles approach each 
other. Thus the above division by the factor Y[i i<ji'^i ^ 
Zj) will never produce any poles in the resulting function. 

Therefore to classify FQH phases of electrons, we can 
restrict our attention mainly to symmetric polynomials 
^{{z{}), where $ is invariant under the interchange of 
zf and Zj for any i and j, but not necessarily invariant 
under the interchange of z/ and z'J if / ^ J. In this paper 
we will often refer to / as a layer index. In the follow- 
ing, <I> will always refer to such a symmetric multilayer 
polynomial. 

We will introduce data, such as n, m, and Ss, to 
characterize and classify bosonic FQH states (ie sym- 
metric polynomials) $({z/}). From the above discus- 
sion, we see that the same set of data also characterizes 
fcrmionic FQH states whose wave functions are given by 

^a„ti-symm({4}) = '^{{4})Ul;i<j{zi " Zj)- 

III. IDEAL HAMILTONIANS AND IDEAL 
WAVE FUNCTIONS 

Progress in understanding the various possible topolog- 
ical phases in the FQH effect has occurred largely because 
of the discovery of various kinds of model wave functions 
and the ideal Hamiltonians that yield these wavefunc- 
tions as their ground states. In the pattern of zeros ap- 
proach, we classify all of the possible ideal wave functions 
for FQH phases. In this section we will explain what we 
mean by "ideal." 

For each topological phase in the FQH system, we want 
to have a certain representative wave function that cap- 
tures the topological properties of the phase. The pro- 
totypical example of such an ideal wavefunction is the 
Laughlin wave function at filling fraction v = 1/m: 

^l/m{{z^}) = Uiz. - Zj)^. (4) 
i<j 

At every particle in there are m zeros, and there 

are no off-particle zeros. The exact ground state of the 
microscopic Hamiltonian with the Coloumb interaction 
does not have this simple property. Nevertheless, the 
ideal Laughlin wavefunction captures all of the essential 
topological properties of these phases. 

The Laughlin wavefunctions arc also ground states of 
certain ideal Hamiltonians that contain interaction po- 
tentials that are constructed only from J-functions and 
their derivatives. For example, if the interaction poten- 
tial is of the form 

Vi{zi,Z2) = 5{zi - Z2), (5) 
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then the following wavefunction is an exact ground state, 
with vanishing total potential energy: 



*i/2 = Y[{z^ - Zjf 



(6) 



On the other hand, if the interaction potential between 
two electrons is of the form 



V2{zi,Z2) = VqS{zi - Z2) + V2d%6{zi - Z2)dl 



(7) 



with > and V2 > 0, then the zero-energy ground 
state will be 



*l/4 = WiZi - Zjf. 



(8) 



More complicated ground states can be obtained 
through more complicated interactions. For example, 
consider the following three-body interaction between 
electrons: 



Vpi =S{vo5{zi - Z2)S{Z2 - Z3) 

- vi6{zi - Z2)d;,*5{z2 - Z3)d;,s), 



(9) 



where S is the total symmctrization operator between 
Zi, Z2, and Z3 . An exact zero-energy ground state of this 
interaction is the Pfaffian wavefunction a.t u = 1: 



$Pf = A 



1 



Zi Z-i 



i<j 



(10) 



where A is the total antisymmetrization operator be- 
tween Zl, - ■ ■ , Zj^. 

Since the Pfaffian and Laughlin wavefunctions are ex- 
act zero-energy ground states of interaction potentials 
that arc constructed only from ^-functions and their 
derivatives, these wavefunctions can be characterized 
completely by their pattern of zeros, i.e. by the order 
of the zeros in the wavefunction as different numbers of 
particles are taken together. These wave functions have 
certain ideal properties; for example, the zeros that are 
bound to a particle lie precisely at the location of that 
particle, not slightly away from it. In this paper, we will 
classify ideal FQH wavefunctions, which are wave func- 
tions that are exact zero energy ground states of such 
ideal Hamiltonians. 

The existence of these ideal Hamiltonians is crucial. 
We can write down any arbitrary complex function of TV 
variables, but we cannot know whether it corresponds to 
a valid topological phase of matter in the limit N ^ 00 
unless we also know that it is the ground state of some 
local, gapped Hamiltonian. If this is the case, then it 
is conceivable that there exists some physical situation 
in which the low energy effective interactions between 
electrons yields a phase of matter that is in the same 
universality class as the ideal wave functions that we 
constructed. Unfortunately, judging whether a many- 
particle interacting Hamiltonian is gapped in the ther- 
modynamic limit is in general an intractable problem. 



So here we limit ourselves to classifying ideal wave func- 
tions, which at least we believe can be realized as ground 
states of local Hamiltonians. Whether the corresponding 
Hamiltonian can be gapped in the thermodynamic limit 
is a question that we must attempt to answer along the 
way. 

The ideal quasihole wavefunctions for these states are 
also zero energy cigenstates of their corresponding ideal 
Hamiltonians. The true ground state is distinguished by 
being the unique translationally invariant state with low- 
est angular momentum. Thus we can also classify topo- 
logically distinct quasiholcs by their pattern of zeros, i.e. 
by the order of the zeros in the quasihole wavefunction as 
different numbers of electrons are taken to the quasihole. 
While the quasihole wavefunctions are also zero energy 
states, the quasiparticles will be gapped and will have the 
same topological properties as the quasiholes. Therefore 
in the following, when we discuss the pattern of zeros of 
the quasiparticles, we are referring to the pattern of zeros 
of the quasihole wave functions. 

Not all FQH phases have such ideal wave functions. 
The hierarchy states and the composite fermion states, 
for example, do not have ideal wave functions. While 
these phases do have their own model wave functions, 
they are not ideal in the sense that they cannot be di- 
rectly described by their pattern of zeros or written as 
a correlation function of conformal primary operators in 
a CFT. Therefore the pattern-of-zeros construction does 
not directly classify these phases. We will however dis- 
cuss how they are related to the ideal multilayer FQH 
wave functions in section IX A. 



IV. PATTERN OF ZEROS 
CHARACTERIZATION 

The spirit of the pattern-of-zeros approach is to con- 
sider bringing together a/ particles of type /, for / = 
1, • • • ,Nf, and asking how $ goes to zero under such a 
procedure. The order of the zero will be denoted S3, 
where a = {ai, - ■ ■ ,aNf)- In the following we will more 
precisely define Ss and discuss some different yet equiv- 
alent ways of characterizing the pattern of zeros. This 
discussion is a straightforward generalization of the dis- 
cussion in the single-layer case.^° 



A. Sa characterization 

Consider a set of a/ coordinates of each type /, and 
set a = (fli, • • • , ftAT^). Define Ss as the minimal power 

of (nfii nr=i zl) in the polynomial This means that 
if we set 

z'i = X^'i + z^^^ i = l,---,ai, V/, 

z(^ = hi^^^, Y.^l = ^, (11) 
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and we take A — > 0, then 

$ ~ X'^Pim, • • • }) + 0{X'-+'), (12) 

where {^ar+i' • • • }) is a polynomial in {^/} 

and the remaining coordinates z*-"-* and {z^^_^_j}. We refer 
to z^") as the coordinate of an a-cluster. We assume that 
Sg is independent of the choice of z^") , which must be the 
case for translationally invariant wavcfunctions. Wc also 
assume that Sg is independent of the choice of {^/} and 
that different polynomials -P({^/}, z^°-\ {2^7+1; ' ' ' }) ob- 
tained from different c;lioices of ^/ arc linearly dependent. 
This is the assumption of unique fusion . 

We can immediately deduce some basic properties of 
Sg. Since $ has no poles, it is clear that Sg > 0. Since 
$ must be single- valued under rotating A in the complex 
plane by an angle 27r, Sg must be an integer. Let Sgj 
be the minimal power of z{: that is, (ei).j = Sjj. A 
translationally invariant $ will have Sgj = 0, otherwise 
it will vanish everywhere. 

Thus, for a translationally invariant polynomial, Sg is 
a nonnegative integer that characterizes the order of zero 
that results when the size of an a-cluster goes to zero. 

B. Derived Polynomials and the D-j 
chciracterization 

In the previous section, we introduced the derived 
polynomials -P({^/}, {2a/-i-i> ' ' ' })• a conse- 

quence of the unique fusion condition, these polynomi- 
als are actually independent of {^/}. We may consider 
more general derived polynomials by bringing together 
other sets of coordinates in P to obtain P{z^'^\ z^''\ ■ ■ ■ )■ 
Then we may consider bringing together an a-cluster and 
a 6-cluster: 

^{Z ,Z , >\^(S)_^(b)^^(a + b) 

+ 0((^(") (13) 

Thus, D-r characterizes the order of the zeros in the de- 

' ab 

rived polynomials as a cluster of a electrons arc brought 
close to a cluster of b electrons. The unique-fusion condi- 
tion assumes that the derived polynomials obtained from 
different ways of fusion are always linearly dependent. 

The fact that $ is a single-valued, symmetric polyno- 
mial implies 

Dsb = Dss&^, Daa=eyen, D,^>0. (14) 

Wc can deduce a relation between D-g and Sg as fol- 
lows. The order of the zero obtained by creating an 
{a + 6)-cluster is -S-^g. One way of creating such a cluster 

is by first creating an a-cluster, then creating a 6-cluster, 
and finally bringing together the two clusters to create 
an (a -t- 6)-cluster. The order of zero in this case will be 



Sg + + D-g. Thus D-j can be obtained from {Sg} 
through the formula 



Since Sg^ = 0, where recall (e/)j = Sjj is the unit vector 
in the / direction, we also have 

Sg+ei =83 + Dg^gj . (16) 

From this recursion relation and from the fact that Sg^ = 
0, we may obtain Sg from the sequence -D-g. Therefore 
we may equivalently label the pattern-of-zeros data using 
{Sg} or {D.j;}. 

C. Characterization by sequence of highest 
occupied orbitals 

The integer Sg has the following meaning. A polyno- 
mial with a I particles of 7**^ kind {ie aj particles in the 
layer) has a total order of Sg. In other words the 
total angular momentum of the quantum Hall droplet is 
Sg if the droplet has particles in the 7*^ layer. If we 
remove one particle from the 7*^ layer, then the total 
angular momentum of the quantum Hall droplet will be 
reduced to Sg-gj . Thus we can interpret 

Ig = Sg Sg—gj (17) 

as the angular momentum of the highest occupied orbital 
in the 7**^ layer for a quantum Hall droplet with a"^ par- 
ticles in the J*^ layer. The iV/-dimensional sequence of 

vectors Ig = {lg,--,lg^) will be called the sequence of 
highest occupied orbitals (HOO). 

We see that ll makes sense only when > 0. Wc 
will set Ig = when it does not make sense. From (17), 
we also see that there is one-to-one correspondence be- 
tween the sequence Sg and Ig. Thus wc can also use Ig 
to characterize the pattern of zeros in the wave function. 

D. Relation to Angular Momentum on the Sphere 

A FQH wave function ^{{z(}) defined on a sphere 
forms a representation of SU{2). In such a case, z rep- 
resents the stereographic projection onto the plane of a 
point on the sphere. A single particle in the lowest Lan- 
dau level can fill any of the 7V$ -|- 1 orbitals, so the rep- 
resentation of SU{2) formed in this case is the one with 
angular momentum J = iV$/2. The SU(2) Lie algebra is 
generated by 

= zd, - J, L- =d,, L+ = -z^d^ + 2Jz. (18) 

In the multilayer case, then, the angular momentum of 
a particle of type 7 will be J/ = where is the 

total number of flux quanta through the sphere seen by 
the particles in the 7"^ layer. Note that here we allow the 
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numbers of flux quanta in different layers to be different. 
The total angular momentum of an a cluster in the z 
direction will be the eigenvalue of the operator 



(19) 



i=l 



The operator ^^lY^lLi^idz' counts the total power 
of a polynomial. Since the minimum total power of 
11/ nr=i the minimum total angular momentum 

of an a-cluster is given by 5's — X]/ cliJi. This means that 
the a-cluster carries an angular momentum of 



Js = a ■ J - Ss = -a ■ Nis, - Ss, 



(20) 



where J = (Ji, Jjv^) and = {N^, 0- We will 
use this relation later to construct ideal Hamiltonians and 

to place conditions on the pattern of zeros for when they 
can correspond to rotationally invariant wave functions 
on the sphere. 



in order to diminish the strength of the on-particle 
zeros. By definition, the phase change of / under the 
above procedure is 27rZ)-j^-. Therefore, the condition 
that the wavefunction have no poles leads directly to the 
following concavity condition on the integers D^p 



a,b+c — ab ' 



(21) 



In cases where all of the zeros arc located on the par- 
ticles and there are no off-particle zeros, the above in- 
equality is saturated. This occurs in the Laughlin states 
<[> = ni<j(-^i ^ ^^"^ their multilayer Abelian gener- 

alizations, the Halperin states 

* = n - ^jf" n - ^if"- (22) 

I-.i<j I<J;i,j 



In the following we will rewrite the concave condition 



as 



A3(a,6,c)>0, 



(23) 



V. CONSISTENCY CONDITIONS 

For the pattern of zeros to describe a valid FQH wave- 
function, it must satisiy certain consistency conditions. 
We already encountered several such conditions above. 
For instance, we found that Ss is a nonnegative integer, 
^ab = ^a+b - Sa- S^> 0, and Dss = S2S - is even. 
In the following we find additional conditions that the 
pattern of zeros must satisfy. 

A. Concave condition 

One of the most important conditions on the wave 
function is simply the condition that the wave function 
have no poles. This condition is remarkably restrictive 
on the allowed pattern-of-zeros sequences. Consider a de- 
rived polynomial P{z^'^\ z^''\- ■ ■ ) and fix all coordinates 
but z^°-\ thus viewing it as a complex function /(z*^"^). 
/(z^")) has zeros at isolated points, but no poles any- 
where. Some of the zeros are located at z^''\ z^^, etc . 
Those zeros arc called on-particle zeros. The rest of the 
zeros are called off-particle zeros. 

If we imagine taking z^"-^ around z^''^ without enclos- 
ing any off-particle zeros, then / will pick up a phase 
27r£'-g. Similarly, if we take 2;^"^ around without en- 
closing any off-particle zeros, then / will pick up a phase 
2ttDss- Now consider taking Under such a 

process, some nearby off-particle zeros will also be taken 
to z^^ . Therefore, if we take z^"^^ around a contour that 
encloses both z^''^ and z^^ in the limit that z^''^ z^^ , 
the complex function / must change by a phase that is 
greater than or equal to 2-k{D^^ + Dg^). The phase can 
never be less than this amount because that would re- 
quire the existence of off-particle poles that get taken to 



A3(a, b, c) ^D,^^, D,^ + Ds^ 

+ Ss + S-, + Ss. (24) 

B. Cluster Condition 

The cluster condition is a way to associate some kind 
of grading to the polynomials that is physically mean- 
ingful. Let {n/} for J = 1, • • • , A^/ be a set of vectors 
that generate an A^/-dimensional lattice, where Nj is, as 
before, the number of flavors of particles (or the number 
of layers). The cluster condition states that the concave 
condition is saturated, i.e. 

DaUc = ^Sb + Dac, (25) 

if either a, 6, or c lie on the lattice generated by {n/}. 
That is, if either a, b, or c can be written as a linear 
combination with integer coefficients of the vectors {n/}. 
This means that a derived polynomial containing a k = 
X]/ kifii cluster is non-zero unless z^''^ coincides with 
the coordinates of another cluster; viewed as a function 
of the single variable z*^'^^ it has no off-particle zeros. A 
consequence of this is that if all of the particles are fused 
to form n-clusters, then the resulting derived polynomial 
has the Laughlin-Halperin form (see (22)) and there are 
no off-particlc> zeros. 

The siugle-layer Read-Rezayi Zn parafermion wave 
functions satisfy an n-clustcr condition and they are ex- 
act ground states of Hamiltonians with n + 1-body in- 
teractions. For a fixed filling fraction, as n increases, 
the number of topologically distinct quasiparticles, the 
ground state degeneracy on higher genus surfaces and the 
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complexity of interactions necessary to realize the state 
all increase. This suggests that the energy gap typically 
decreases with increasing n. Wave functions that do not 
obey a cluster condition can be thought of as having in- 
finite n and are not expected to correspond to gapped 
phases. This intuition also comes from the CFT approach 
to FQH wave functions; infinite n corresponds to an ir- 
rational conformal field theory, which does not yield a 
finite number of quasiparticles and a finite ground state 
degeneracy on the torus. In the multilayer case, we may 
use the volume of the unit cell spanned by {n/} as one 
way to measure the complexity of a given FQH state. 

The cluster condition is extremely powerful and simpli- 
fying because it allows us to determine the entire pattern- 
of-zeros sequence from knowledge of a "small" number of 
them. To see how this works, first observe using (25) 
that 

Dnj,S = XI ^I^nj,ei = '^JI^^I' (26) 

where we have defined the matrix mji = Z?fj,,(Tj. So 
for any vector k = J2i kinj, where kj is an integer and 
X^/ ki{ni)j > 0, we have: 



(27) 



ij 



The above equations imply 
Dfiunj = ^njADfii,eA = ^njAmiA = {nmF)ji 

A A 

= ^niADg^^fij = niAmjA = {mn^)ji, 

A I 

(28) 

where we have also defined the matrix ri/j = {ni)j. 

In terms of the sequence {Ss}, this implies that for 
k = ^jkini, where ki is an integer, 

= Sa + ^ kiSfij + ^ kimijaj 



IJ 



+ ^ ^{nm^)ji{kikj - 5ijki). 



(29) 



IJ 



Therefore, all of the integers Ss are specified by the points 

a within the unit cell spanned by {n/}. 
In terms of the HOO squence Is, we have: 
/^^— Q <? ^ 

S+k a+k '-'a+k-ei 



= Ss - Ss-ei + ^ kAmAI 
A 

= Is + ^^ArriAi- 



C. Equal Area Layers 

The density profiles of the single-particle states in the 
lowest Landau level, 2;"*e~l^l /^'b, are in the shape of a 
ring, with a peak at a radius = V^mls, where Is is 
the magnetic length. Such a wave function has an an- 
gular momentum m. When many of these orbitals are 
filled by particles, the total wavefunction will describe a 
uniform, rotationally symmetric state that goes to zero 
at a radius Vmax = V^JT^maxls, where mmax refers to the 
filled orbital with maximum angular momentum. There- 
fore, a given quantum Hall wavefunction will describe a 
QH droplet of area Aimimaxl'B, and rrimax is given by the 
maximum power of zi (or Zi for any other fixed i) . rrimax 
is also equal to the number of flux quanta, . 

An important constraint on the multilayer quantum 
Hall wavefunctions is that they must describe systems 
in which each layer occupies the same area, up to small 
corrections. The requirement that each layer occupies 
exactly the same area amounts to the requirement that 
each layer has exactly the same number of flux quanta, 
N^, = = N^. However, it is reasonable to include 
states in which different layers occupy equal areas only 
up to 0{N°) corrections, where is the number of elec- 
trons. 

Such a requirement of approximately equal area layers 
is summarized in the following equation: 
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1. 



(32) 



We wish to see how this condition translates into a condi- 
tion on the pattern of zeros. The conditions are slightly 
different depending on whether we ultimately want to 
characterize gapped FQH phases of fermions or bosons. 
If we are interested in fermionic phases, we require that 
= Ui-i<j{zi - z'jMW}) be a valid FQH wave 
function of fermions, which does not reqiure that ^{{zf }) 
be a valid FQH wave function of bosons. In what follows 
we will explicitly analyze the bosonic case, where we re- 
quire $({2:/}) to be a valid multilayer FQH wave function 
of bosons. 

is equal to the maximal power of zf; for the boson 
wavefunction this is given by 



(33) 



where recall Nj is the number of particles of type /. Us- 
ing the cluster condition, we find 



(30) 



Ni = Y, Ni{n-^m)ij + Sftj - S^j 



mjj, (34) 



Finally, note that since Ds^s is an even integer, we have: 
even = Dfijnj = '^njiDfij,ei = ^njimji 



= {nm^)jj. 



(31) 



where we have set N = '^jkinj. Requiring (32), we 
obtain the following condition on the pattern of zeros: 



(35) 
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This can be seen most easily by ignoring the 0{N^) terms 
in (34), taking N^/N^ 1, and inverting the result to 
obtain Nj ~ N^'^j(m~^n)ij, which must be nonneg- 
ative. From this analysis, we learn that if n~^m is not 
invertible, then the pattern of zeros cannot fix the ra- 
tio of particles Nj/Ni in the different layers. Therefore 
the corresponding FQH state has a gapless mode corre- 
sponding to the relative density fluctuations between the 
different layers. 

As a simple example of this analysis, consider the 
(1, 1, 1) Halperin bilayer state, which is known to have 

a gapless density mode and for which n~^m — ^| | 

which is not invertible. A macroscopic number of par- 
ticles can freely go from one layer to the other with- 
out changing the area of the quantum Hall droplets, sig- 
nalling the existence of a gapless relative density mode. 
Inverting (34) yields 

- + Sn,-ej + mil) = Nj, (36) 



from which we can read off the filling fraction in each 
layer: 



1^1 



E 



(to ^n)ij. 



(37) 



The total filling fraction is the sum of the filling fraction 
of each layer: u = i/j. 

For fermions, (33) is modified to 



Ni = S^-S^_,,+Ni 



N-ei 



(38) 



due to the extra factor n/j<j(-^/ ^ ^ii^i})- Note 

that {Ss\ still describes the pattern of zeros of the sym- 
metric polynomial i>. The result for fermions is therefore 



VI = ^(I + n-^m)j} > 0, 



(39) 



where <S is of order 1 in the large A^e limit (see eqn. (36) 
)• 

More precisely, completely filling the sphere means 

that the quantum Hall state is rotationally invariant with 
zero total angular momentum. Using (20), we find that, 
for a bosonic FQH state characterized by Ss with A/ par- 
ticles and flux quanta in the I^^ layer, the maximum 
total angular momentum is given by 



(41) 



where N = (Ai, Ajv^) and A* = (Ai,...,A^O- For 
a fermionic FQH state characterized by Sg with A/ par- 
ticles and A| flux quanta in the I*"^ layer, the maximum 
total angular momentum is given by 



J^ = -N.N^-S^ 2 



In the above, A$ must satisfy 



A| > Sf} Sj^^g^ 



IY,Ni{Ni-1). 



(42) 



for bosons 



A| > 5^ - +N1-I for fermions 



(43) 



in order for the wave function to fit into each layer. 
Completely filling the sphere requires that 



for bosons 
A/ — 1 for fermions 



(44) 



and Jjq = 0. We see that A and A$ must satisfy 



A- A$ = 



25*^^ for bosons 

2Sj^ + J2i {Ni-1) for fermions 



(45) 



where I is the A/ x Nf identity matrix. If (I -|- n~^m) is 
not invertible in the fermionic case, then there are gapless 
relative density modes, which is why the filling fraction 
in each layer becomes undefined. 



which implies (for both bosons and fermions) 

Y^Ni{S^-S^_.^) = 2S^. (46) 



D. Shift and Rotational Invariance on the Sphere 

Consider a multilayer quantum state with Aj particles 
in the 7**^ layer. We want to put the quantum state 
on a sphere with A^ flux quanta in the layer. We 
would like to know for which set of A^ can the quantum 
Hall state completely fill the sphere? Naively, one may 
expect A^ and A/ are related by the filling fraction in 
each layer N^/i/i = N^. However the precise relation 
between the number of flux quanta and the number of 
electrons includes a shift, 



j/-i^z//A| =!/-^Ae-5, 



(40) 



If a given A docs not satisfy (46), then the corresponding 
quantum Hall state (with A/ particles on the I^^ layer) 
cannot completely fill the sphere. For A that satisfies 
(46), the corresponding quantum Hall state can com- 
pletely fill the sphere and has zero total angular momen- 
tum if A$ is given by (44). (46) can generally be satisfied 
only if A lies on the lattice spanned by {n/}. 

We would like to remark that it is easy to have different 
numbers of flux quanta on different layers in numerical 
calculations. The pattern of (A, A$) where the quan- 
tum Hall state has zero total angular momentum on the 
sphere can be used as a fingerprint to identify different 
quantum Hall states through numerical calculations (for 
examples, see Tables IV and V). 
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E. Additional constraints: A3 = even 

The analysis of the single-layer case in Ref. 10 has sug- 
gested an additional condition: 

A3 (a, 6, c) = even. (47) 

There, it was found that allowing A3 (a, 6, c) = odd allows 
for certain pattern-of-zeros sequences that either do not 
correspond to single-valued wavefunctions (such as the 
square root of the PfafRan) or could not correspond to 
translationally invariant wavefunctions. It was suggested 
that one way to rule out such possibilities is to impose 
(47). How should this condition be generalized to the 
multilayer situation? 

One natural generalization is to impose A3 (a, b,c) = 
even for all o, b, and c. However, we find that this con- 
dition is too restrictive. It rules out certain known FQH 
wavefunctions, such as the sw(3)2/w(l)^ non-Abelian 
spin singlet states. ^^'^^ The need to relax this condition 
while still having it remain compatible with the single- 
layer situation suggests that we should impose (47) only 
for choices of a, b, c that are collinear through the origin. 

While it was found that allowing A3(a,5, c) = odd 
allows for pattern-of-zeros sequences that do not seem 
to correspond to valid translationally invariant, single- 
valued wavefunctions. there arc known cases of CFTs 
with A3 (a, b, c) = odd that do seem to yield translation- 
ally invariant, single-valued wavefunctions. One such ex- 
ample is the so-called Gaffnian wavefunction, which has 
A3 (1,1,1) = odd and which can be constructed using 
the minimal model CFT M{5, 3)}^ This CFT however is 
non-unitary. It has been suggested that FQH wavefunc- 
tions constructed using non-unitary CFTs correspond to 
gapless phases;^*^ whether this is always necessarily the 
case is currently an important open question in FQH the- 
ory 

These considerations suggest that in order to restrict 
ourselves to pattern-of-zeros sequences that have a cor- 
responding unitary CFT, we should impose A3 (a, b,c) = 
even for those a, b, and c that are collinear through 
the origin. In our search for pattern-of-zeros solutions, 
we will impose this condition and analyze the resulting 
states. The precise connection, if any, between this con- 
dition and valid FQH wavefunctions that correspond to 
unitary CFTs remains to be clarified. 



VI. IDEAL HAMILTONIANS 

Given a pattern-of-zeros sequence, it is important to 
be able to construct a local, gapped Hamiltonian whose 
ground state wave function has the given pattern of zeros. 
If this is possible, then we know that the corresponding 
pattern-of-zeros sequence describes a topological phase of 
matter. Whether this particular phase is realized in an 
experiment then depends on the particular types of low 



energy effective interactions between the electrons in the 
unfilled Landau levels. 

We can go about constructing such a Hamiltonian by 
noticing that on a spher(\ the integers Sg are directly 
related to the angular momentum of the a-cluster. For 
an electron system on a sphere with A^^ flux quanta for 
the Ith layer, an electron of type / will carry an angular 
momentum = N^/2. For an a-cluster, the maximum 
angular momentum is therefore a ■ J. However, for a 
polynomial ^{{zj }) described by a pattern of zeros {Sg}, 
the maximum allowed angular momentum of the a cluster 
is only = a ■ J — Ss- The pattern of zeros forbids the 
appearance of angular momentum a-J-^Ss + ^ja-J — 
Ss + 2, - ■ ■ , a • J for any a-clusters in c&dz-}). 

Such a condition can be easily enforced by writing the 
Hamiltonian as a sum of projection operators, -Pj^^ Let 
Pg"-* be a projection operator that acts on the a-cluster 

Hilbert space. Pg"^ projects onto the subspace of a- 
clusters that have total angular momenta greater than 
a - J — S. Now consider the Hamiltonian 

a a-clusters 

where X^s-ciustcrs sums over all of the a-clusters for a fixed 
o. The wave function described by {S^} will clearly be 
a zero-energy ground state of the above H^Ss}- In many 
cases, there is only one unique ground state wave func- 
tion with minimal total angular momentum, but in gen- 
eral there can be many independent polynomials with the 
same pattern of zeros. In such a situation, the Hamilto- 
nian would need to be modified further to select for a 
particular polynomial with the given pattern of zeros. 

In order for the above Hamiltonian to be local, J^s 
must be limited to a small, finite number of a-clusters. 
But as a result, we cannot guarantee in general that the 
ground state wave functions will all be described by the 
sequence {Sg} for every a, or even that they will obey the 
cluster condition. In many of the known cases, such as 
the Laughlin, Moore-Read, and Read-Rezayi wave func- 
tions, the sum over a-clusters can indeed be terminated 
after the first few clusters while still yielding a unique 
zero energy ground state wave function with minimal an- 
gular momentum which is described by {S's}. 

The above construction for H^Sg} should therefore 
be viewed as a starting point for constructing an ideal 
Hamiltonian that is local, gapped, and whose unique 
minimal angular momentum ground state wave function 
is described by {Sg}. In some of the simplest cases, we 
know that this construction suffices. 



VII. SUMMARY: PATTERN-OF-ZEROS DATA 
AND CONDITIONS 

We have found that the polynomials $({z/}), / = 
1, • • • ,Nf, that may correspond to stable FQH states are 
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described by the following data: 

n, m, {Ss}, (49) 

where n and m are Nf x Nf matrices with integer entries 
that satisfy 

mij > 0, nij > 0, det n^O 

mri^ = nm^ , {mn^)ii = even. (50) 

The above implies that n~^m is a symmetric matrix. 

Furthermore, for the pattern of zeros to fix the relative 
densities of particles in each layer, we have 

{n^^m)ij is invertiblc (for bosons) 
{n~^ni)ij + Sjj is invertiblc (for fcrmions). (51) 

Otherwise, there are gapless relative density fluctuations. 
We also have 

ui > 0, 

I, = I I2Am-^n)ij ^ for bosons 

\ X^j(I + n~^m)jj for fermions ^ ^ 

The total filling fraction is u = vi. 

Note that the {Ss} need to be specified only for values 
of a that are contained in the unit cell spanned by {n/}, 
where nj corresponds to the 7*^ row of the matrix n: 

{ni)j = nij. 

Using the cluster condition, we can determine Ss for 
all a simply from m, the fact that the Sg^ = 0, and from 
the values of 5g for which 6 lies in the unit cell spanned 
by {n/}: 

Ss_^_k = Ss + ^ kiSfij + ^ kimijaj 
I ij 

+ ^ y^\nm^)ij{kikj - dijki), (53) 
IJ 

where k = kinj. Ss must satisfy: 

A.2{a, a) = even, 
A2(a,&)>0, A3{a,b,c)>0, (54) 

where 

A2(a,6)^5,^g-55-5g, 

A3 (a, b, c) =S^^^^^ - S^^^ - Ss^o - -S^b+e- 

+ Ss + S-^ + Sn. (55) 

Finally we impose 

A3 (a, b,c) = even (56) 
for those a, 6, and c that are coUinear through the origin. 



A. Product of Symmetric Polynomials and 
Primitive Solutions 

Consider two polynomials of r7-cIuKtcr form, $ and <!>', 
and consider their product: $ — The pattern of 

zeros of $ is the sum of the pattern of zeros of $ and $': 

a/3 af3 ' 5/3 ^ ' 

Similarly, the data in terms of m and Ss are also additive. 
Note that all of the conditions that we impose on the 
pattern-of-zeros are linear. Thus the pattern of zeros of 
$ is also valid. However, notice that the condition for 
filling the sphere is not linear. There may be two FQH 
wavefunctions $ and that can fill the sphere but whose 
product $ cannot fill the sphere. 

Thus, we can divide the pattern-of-zeros solutions into 
primitive and non-primitive solutions. Primitive solu- 
tions are solutions that cannot be written as a sum of 
two other solutions. 



VIII. RELATION TO CONFORMAL FIELD 
THEORY 

The pattern-of-zeros approach is closely related to the 
conformal field theory approach to constructing FQH 
wavefunctions. In the CFT approach, the symmetric 
polynomial $({z/}) that describes a multilayer FQH 
state can be written as a correlation function of a set 
of electron operators Vg-j in a CFT:^^ 

$({4})= lim z'J^^{V{z^)llVe.Az()). (58) 

i,I 

The operators Vg-j are written in the form: 

Ve.j{z) = i;g,{z)e'^^^"1>'^'\ (59) 

where e^^J ^^'-''^•'^^^ is a vertex operator in a U{l)^f 
CFT. It has scaling dimension Mfj/2. is a simple 
current operator; that is, it satisfies the following fusion 
relation: 

V-sV-s = i^s+b- (60) 

This Abelian fusion rule is the CFT version of the unique- 
fusion condition. The cluster condition implies that ip 
satisfies 

^n. ~ 1, (61) 

where (n/)j = nij. An a-cluster of electrons will be 
described by the operator 

Thus we see that the cluster condition implies that 
an n/ cluster is described by a vertex operator 
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e*E,7K"/.7MjK*K(2)^ If all of the particles are grouped 
into n-clusters, then the corresponding derived polyno- 
mial will be just a correlation function of vertex opera- 
tors in a U{l)^f theory, which will have no off-particle 
zeros and will be of the Laughlin-Halperin form. 
Let us denote the scaling dimension of the operator V3 

as 



ha - ht^ + hf 



(63) 



where h'^g is the scaling dimension of the simple current 
ips and h^'^ is the scaling dimension of the vertex opera- 
tor. Notice that since ■f/'rij ^ 1, the simple current scaling 
dimensions satisfy /iS^- = h^^. The scaling dimension 
of the Gaussian part is given by 



hf = -ai{MM^)i.jaj. 



(64) 



The pattern of zeros are related to the scaling dimensions 
through the relation 



(65) 

This allows us to obtain the scaling dimensions from the 
pattern of zeros. Using the cluster condition, some alge- 
bra shows that MM^ = n~^m and so (64) becomes 



hf- 



(66) 



The scaling dimensions of the simplc-ciirrent part can 
also be determined from the pattern of zeros by using 
the fact that hftj = (ihilJii^ fi^ = q, and applying (65) 
iteratively. This yields: 



II 



= '^niAhg^ + Sfij. 



(67) 



Thus the FQH wave function is of the form 

$({4}) = $,e({-^/})$ga({-^/}), 

where 



(71) 



^ga({4}) = 

JJ (^[ - 4)(""''")" II (4-^/)(""'™)". 



I;i<j 



I<J;i,j 



(72) 



$sc arises from the correlation function of the simple cur- 
rent sector and is the "non-Abelian part" of the wave- 
function. 

In this manner, each pattern-of-zeros solution corre- 
sponds to the current algebra of a rational CFT. The 
connection between the pattern of zeros approach and 
the CFT approach can be thought of in the following way. 
The pattern of zeros describes the essential properties of 
the CFT that yield valid FQH wavefunctions. So in or- 
der to classify ideal FQH wavefunctions, one can bypass 
the CFT altogether and go directly to the heart of the 
matter: characterizing the allowed pattern-of-zeros solu- 
tions. Furthermore, since each pattern-of-zeros solution 
corrresponds to a CFT, the pattern-of-zeros classification 
can be viewed as a classification of the allowed CFTs that 
can be used to construct FQH wave functions. 

In this formulation, the pattern of zeros classifies all 
those ideal FQH wavefunctions that can be formulated as 
a correlation function of conformal primary fields. There 
are also many FQH wavefunctions, such as the hierarchy 
states and the Jain scries, that cannot be written in this 
way. These wavefunctions are outside of the pattern-of- 
zeros classification. However they may be closely related 
to the pattern-of-zeros construction. We comment on 
this connection elsewhere. 



Multiplying both sides by n ^ gives 

h,^ = hZ + = n-( - S,^ (68) 

In a similar manner, one can obtain 



hs = ^ai{n ^)ij{ 



h"/ 



-T -1 ^ 

a n ma 



(69) 



I:i 



which determines h'^g in terms of the pattern-of-zeros. 

Note that the correlation fmiction of the Gaussian part 
is, leaving the background charge implicit, 

= n - 4)^'''''^" n (4-^/)^^^"'" 

I;i<j l,J;i,j 

(70) 



A. Alternate Labeling 

Using (69), we can derive a formula for Ss in terms of 
h%'^, n, and m: 



Sa = hs- ^ aihgj . 



(73) 



Thus there is a one-to-one correspondence between the 
simple-current scaling dimensions and the sequence {Ss}- 
This means that there is yet another way to label the pat- 
tern of zeros. Earlier, wo found that one convenient label- 
ing of the pattern of zeros is with an Nf x Nf matrix n, 
an Nj X Nf matrix m, and the value of the non-negative 
integers Sa for a lying inside the unit cell spanned by 
the rows of n. An alternative, equivalent labeling of the 
pattern of zeros is by specifying the following data: 



n, m. 



(74) 



for a lying inside the unit cell spanned by n. This labeling 

is convenient because it makes close contact with the cor- 
responding CFT description; if h^S ^ for some a, then 
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the CFT has a non-trivial simple-current structure and 
therefore generally also has some form of non-Abelian 
statistics. 



B. Relevant CFTs for multilayer FQH states: 
gk/u{iy parafermions 

In the single-layer case, many of the pattern-of-zeros 
solutions were found to be closely related to the Zk 
parafermion CFTs of Zamalodchikov and Fateev.^^ What 
are the relevant CFTs to expect in the multilayer case? 
The answer is that some of the corresponding CFTs in 
the multilayer case will be closely related to the gk/u{iy 
parafermion CFTs that were constructed by Gepner, 
where gk is a simple afEne Lie algebra at level k and 
r is the rank of the Lie algebra gP'^ The case g = su{2) is 
equivalent to the Zk parafermion CFTs of Zamalodchikov 
and Fateev. 

The simple-current algebra of the gk/u{lY 
parafermion CFT has the following structure. For 
every element a of the r-dimensional root lattice 
of g, associate a simple current operator ipg. The 
simple-currents will have the fusion rules 



(75) 



Furthermore, ips = i>0 H a — P is an element of k times 
the long root lattice of g. The scaling dimension of ips is 
given by 



(76) 



where n{a) is an integer equal to the minimum number 
of roots from which a is composed. The inner product 
= (d, a) is defined with respect to the quadratic form 
matrix of g. 

If wo arc considering quantum Hall states with Nf 
layers, then we would expect to see the appearance of 
these parafermion CFTs with rank r > Nf. Therefore 
in the bilayer case, one class of states that we expect 
to see should be related to gfc/u(l)^ parafermion CFTS 
where 17 is a simple Lie algebra of rank 2. There are 
only three simple Lie algebras of rank 2: su{2>), so(5), 
and G2. Of these, only su(3) is simply laced, so the long 
root lattice is the same as the root lattice. This means 
that for the pattern-of-zcros solutions that correspond to 

, we expect to sec solutions that correspond 



k^ 

to sw(3)fc/u(l)^ CFTs. In Appendix B we will describe 
the sw(3)2/m(1)^ CFT in more detail. 

The parafermion CFTs for g = ,so(5) and g = G2, on 
the other hand, are more complicated because the long 
root lattice is different from the root lattice. For exam- 
ple, so(5)fc/u(l)^ CFTs will generically be relevant for 

2k) (G2)fe/w(l)^ CFTs will generically 

'3k 
3kj 



n 



be relevant for n 



IX. EXAMPLES OF PATTERN OF ZEROS 
SOLUTIONS 

In this section we examine explicitly several simple bi- 
layer pattern of zeros solutions. We first fix the cluster 
structure n to have a simple form. Then we try to find 
all the solutions {83} that satisfy the conditions listed in 
section VII. 

The simplest non-Abelian states can be obtained from 

'1 1' 



the simplest cluster structures n = 



2 



and n = 



2 
2 



These are the simplest non-Abelian generaliza- 



tions of the bilayer Abelian Laughlin-Halperin states. 

Note that by definition the ordering of the rows in n is 
arbitrary; we choose it so that nn ^ 0. Interchanging the 
layers yields the same physical system but corresponds 
to interchanging nn with 7122 and ni2 with 71,21, so two 
matrices n and n' that are related by such an interchange 
are regarded as equivalent. 

To list the solutions for those simplest cases, we may 
use some known CFTs to construct the appropriate 
simple-current algebra that corresponds to the pattern- 
of-zeros solutions. Using this known CFT, we can then 
write the wavefunction explicitly. However, the wave 
function that we write down may not be unique in some 
cases; there may be several independent polynomials that 
have the same pattern of zeros. This corresponds to there 
being several distinct CFTs whose simple-current algebra 
possesses the same pattern of zeros. In the following ex- 
amples, we will make this choice when necessary so that 
we can explicitly write down a wave function with a cer- 
tain pattern of zeros. 



1 
1 



These states are all Abelian and correspond to the 
Halperin (m, n, I) states. The unit cell spanned by the 
nj consists only of the points (1, 0) and (0, 1); By trans- 
lation invariance <S'(i,o) = ^'(0,1) = 0. Thus these states 
are completely characterized by the matrix m and are of 
the form 



(77) 



The m matrix here is exactly the K matrix that describes 
all Abelian FQH states.^ We also have ht^ = for all of 
these states - the CFT has no simple-current part and 
corresponds to a U{1)^ Gaussian CFT. 

Although the single-layer hierarchy states, such as the 
f = 2/5 hierarchy state, do not have ideal single-layer 
wave functions, there are ideal multilayer states that have 
the same topological orders as those single-layer hierar- 
chy states. For example, the topological order in the 
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1/ = 2/5 hierarchy state is described by the ii'-matrix 
(in the symmetric basis). ^'^ Such a topolog- 



K = 



2 3 



ical order can be represented by the ideal bilayer state 
^ ^ ' ^ ' So although the 



with n = 



and m 



Q 1 J y2 3^ 

pattern-of-zcros construction does not directly classify 
those single-layer hierarchy states, their topological or- 
ders can still be described by the pattern-of-zeros ap- 
proach. 



1 1 
2 



This choice of n requires that the electron operators in 
the CFT must take the form 



Vei = i^e'^J , 
Ve2 = i^e'^J^^-"^-', 



(78) 



where ip"^ = 1 and ip has scaling dimension that is integer 
or half-integer. This latter fact can bo obtained from the 
condition A3((l, 0), (1, 0), (1, 0)) = even. One general 
way of constructing such a simple-current operator is by 
expressing it as a set of Majorana fermions from several 
copies of the Ising CFT: 



(79) 



where ^Z'*-*-' is the Majorana fermion from the ith Ising 



CFT. Such an operator has scaling dimension h'^ 
and gives rise to the following FQH wavefunction: 



H{Zi,W,}) = Pi 



X $ 



ga- 



= a/2 



(80) 



Xi represents the coordinates in both layers: 



Xi — 



Zi l<i<Ni 

Wi-N, Ni<i<N2 



(81) 



^ga is defined in (72). Note that the simple-current al- 
gebra in this case implies g) = ^fo !) = -'-• That is, 
these states can all also be viewed as satisfying the clus- 
^2 0\ 

, but with a different choice 



ter condition for n ■ 



2 



m. For every pattern-of-zeros solution found here, there 
is an equivalent one for n = ^ 



2 



Later, we will list the solutions with the n 



2 
2 



cluster structure. Some of those solutions actually have 
il 2)- 



1 1 
1 



Here, the electron operators must be of the form 



V,2 



(82) 



The fact that Hi = (1,1) however also forces ip = 1. Thus 
in fact the corresponding CFTs do not have a simple- 
current part; they all correspond to a J7(l)^ Guassian 
CFT. All of these states are therefore Abelian and corre- 
spond to the Halperin (m, n. I) states. Thus, all pattern- 

Q actually also satisfy 

''1 0^ 

Using the cluster condition and the fact that S'(i,o) = 
^{0,1) = 0, it is easy to see that the pattern of zeros for 
these states arc completely characterized by the matrix 
m and therefore that these states are identical to the 
Halperin states. 



of-zeros solutions where n 



the cluster condition for n = 



D. 



1 
2 



In this case, the electron operators in the CFT must 
take the form 



Vei = e'^-'^i-"^-', 
Ve2 =^e'^-'^"'^-'. 



(83) 



where ^/-^ = 1. From A3((0, 1), (0, 1), (0, 1)) = even, 
it follows that the scaling dimension of ijj is integer or 
half-integer. This means that we can in general write 
it as a product of a Majorana fermion operators from a 
independent copies of the Ising CFT. Thus the pattern- 
of-zeros solution for this choice of n includes only two 
classes of states: the Abelian Halperin states and the 
following version of the Pfaffian states: 



mzi,wi}) = pf 



1 



X 



(84) 



Such a state spontaneously breaks the discrete Z2 sym- 
metry associated with interchanging the two layers. 



E. 



2 
2 



non- Abelian bilayer states 



For this choice of n, the corresponding CFTs have two 

simple current operators, ip(o,i) and V'(i,o); which each 
square to the identity: V'^i 0) ~ V'^o 1) ~ Thus there 
are a total of three distinct primary currents: V'(i,o)) 
■0(0.1); and ip{i,i) = ■0(1,0)0(0.1) • The pattern of zeros can 
be fully specified by specifying the scaling dimensions of 
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these simple current operators, ^^^q-j, ^(oi)' ^^'^ ^'i i)' 
and the matrix m, which specifies the J7(l)^ part of the 
electron operators in the CFT. 

Applying A3 (a, a, a) = even for this choice of n shows 
that the simple-current scaling dimensions are all integer 
or half integer: 2ht'' G Z. One general way of encoding 
these fusion relations and the associated scaling dimen- 
sions is to write the electron operators in the CFT using 
simple-current operators that consist of multiple copies 
of the su(3)2/w(l)^ parafermion CFT. The most general 
choice for the electron operators can be written in the 
form 

Ve;2 = 4r'^ • • • Vg+''"'^e'^^ ^"^^ (85) 

where ipa^ is a simple current from the ath copy of 
the sw(3)2/m(1)^ parafermion CFT and aj and /3i re- 
fer to cither (1,0), (0,1), or (1,1). Some explicit forms 
for such correlators of simple-current operators in the 
S'u(3)2/w(l)^ parafermion CFT were discussed in Ref. 18. 



Computing these correlation functions provides one way 
- not necessarily unique - of constructing a wave function 
with the desired pattern of zeros. 

These pattern-of-zcros solutions can naturally be 
grouped into two distinct classes. In the first class, Ve-i 
and V'e;2 contain only one kind of simple-current, cither 
^"(0,1): ■0(1. o)j or'!/'(i.i), from each copy of the ,s?i(3)2/u(l)^ 
CFT. For example, if '0(i''o) apptiars in Ve-i^ then Ve:2 can- 
not contain or "0(0 !)• such a situation, we can 

think of ipa} as being the Majorana fermion from the 
Ising CFT. This means that these states can be written 
in terms of multiple copies of the Ising CFT . This class 
of multilayer states can be written by choosing 

V;.2 = V'^""'=)---0^°+''"'=)e^^-'^"'^-', (86) 

where ■^^"^ is the Majorana fermion from the ath copy of 
the Ising CFT. The wave function for this class of states 
is therefore: 



(87) 

I 



Xi represents the coordinates in both layers: 
r. = / l<i<Ni 

'''-{wi-N, Ni<i<N2 

This is the simplest generalization of the single-layer 
non-Abclian states to a class of non-Abelian bilayer 
states: the interlayer Pfaffian states. The simplest ver- 
sion of this, with a = 6 = c, is closely related to (and 
identical to for certain choices of m) the spin-charge sep- 
arated non-Abelian spin singlet wavefunction proposed 
in Ref. 24. 

The second class of states cannot be written in terms of 
multiple copies of the Ising CFT - the full su(3)2/u(l)^ 
parafermion CFT is necessary. The first class of states, 
which can be written only using the Ising CFT, have the 
property that their pattern of zeros satisfies A3 (a, b,c) = 
even for all choices of a, b, and c. The second class of 
states, for which the full su(3)2/m(1)^ CFT is necessary, 
satisfies A3 (a, b,c) = odd for certain choices of a, b, and 
c. 

Let us compare the simple-current algebra of the 
su(3)2/u(l)^ parafermion CFT to what one would obtain 
using two copies of the Ising CFT. If we used two copies 
of the Ising CFT, we could have '0(1. 0) = 0''"'^' with scaling 
dimension 1/2, ■0(0, 1) = 0'^^ with scaling dimension 1/2, 



and ■0(1,1) — 0^^V^^^ with scaling dimension 1. This sat- 
isfies A3((l, 0), (0, 1), (1, 1)) = even. On the other hand, 
in the su{3)2/u{l)'^ parafermion CFT, the only difference 
is that ■0(1,1) also has scaling dimension 1/2. Thus in 
this latter theory, two fermions combine to give another 
fermion. This yields A3((l, 0), (0, 1), (1, 1)) = odd. The 
fact that there are valid single- valued translationally in- 
variant FQH wavefunctions that arise from unitary CFTs 
and that have A3 (a, b,c) = odd for certain choices of a, 
b, and c suggests (see Section VE) that we should impose 
A3 (a, b,c) = even not in general but only if a, b, and c 
are collinear through the origin. 

Let us examine the pattern of zeros for a few of the 
simplest examples of these non-Abelian bilayer states. 
There is a fermionic v = 2/3 state with 

'^=(2 2) {'S'(2,o) = 0, S'(i,i) = 0, S'(o,2) = 0} 

u = 2/3 {hll,^ = ^,hll,^ = 0, hll,^ = i}. 

(89) 

This is the pattern of zeros for the interlayer Pfaffian 
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state, which is of the form 



^{{z,,w,})= Pi 
We use the notation 



1 



'^^{2,2,1){{Z^,W,}). (90) 



(91) 



There are also fermionic states at ly = 4/5 and ly = 4/7. 
These have the following pattern of zeros: 



m ■ 



(l 2) {"^(2,0) = 0, 5(1^1) = 0, 5(0,2) = 0} 



u = 4/5 {hll,^ = ^,hll,^ = = i}. 



(92) 



2 3 



[3 2) {"^(2,0) = 0> = 1> '5(0,2) = 0} 

^ = 4/7 {/^a.o) = l^HU = \, HoA) = I}- 



(93) 

The state a.t v = 4/7 is the non-Abelian spin singlet state 
that was proposed in Ref. 17. 

Note once again that the pattern of zeros m and {5s} 
refer to the pattern of zeros of the symmetric polynomial, 
$= * ^ 



X. DISCUSSION OF RESULTS AND 
RELATION TO EXPERIMENT 

In single-layer quantum Hall samples, a quantum Hall 
plateau is seen g,t u = 5/2, but not at u = 1/2. The 
reason is that even though in all of these cases there is 
a single half-filled Landau level, the existence of the two 
filled extra Landau levels modifies the effective interac- 
tions between the electrons in the unfilled level. In the 
J/ = 5/2 case, numerical calculations suggest that the 
these effective interactions are modified in such a way 
that a non-Abelian quantum Hall state may be realized. 

Experiments on multicomponcnt quantum Hall sys- 
tems should be able to probe an even wider variety of 
regimes with distinct effective interactions. For exam- 
ple, for a two-component FQH system, we can study sys- 
tems in which the spin degree of freedom is present, two- 
dimensional electron systems with two quantum wells, 
wide single-layer systems in which the electrons sponta- 
neously form a double-layer system due to Coulomb re- 
pulsion, or systems in which there may be two valleys for 
the free quasiparticle spectrum (such as in graphene or 
SiGe heterostructures), etc. In many of these cases, ex- 
perimentalists can also tune to some extent the degree of 
correlation between the two components. For example, in 



double layer systems, application of a parallel magnetic 
field can tunc the tunneling and correlation between the 
layers. There may also be some degree of tunability in 
the relative densities between the two components in ad- 
dition to being able to probe FQH states with different 
numbers of filled Landau levels. With this greatly in- 
creased amount of variability and tunability in the effec- 
tive interactions between electrons in the unfilled Landau 
levels, it is possible that a non-Abelian state can be re- 
alized in a two-component quantum Hall system. 

Since the pattern of zeros provides a systematic classifi- 
cation and characterization of a wide variety of quantum 
Hall states, it provides us with a general sense of how all 
of the non-Abelian bilayer states are related and which 
ones are simpler than other ones. Just as we know that 
the single-layer Pfaffian quantum Hall state is the sim- 
plest non-Abelian generalization of the Laughlin states, 
we can determine the simplest non-Abelian generaliza- 
tion of the Halperin bilayer states and therefore single 
out some of the possibilities that may be experimentally 
viable. 

In Ref. 25, we have given an overview of some of the 
simplest non-Abelian bilayer states that we find and that 
occur at filling fractions at which experiments on two- 
component FQH systems have already observed incom- 
pressible states. Here we briefly summarize that discus- 
sion and supplement details of the calculations of various 
topological properties of the candidate states. 

Experiments have so far observed FQH plateaus in 
two-component systems at = 2/3, 4/5, 4/7, 4/9, 
6/5, 6/7, 1/4, etc.^^~^^ In some cases, these plateaus 
have been observed in both bilayer and spin-unpolarized 
single-layer systems, while in others, the plateau has only 
been observed in one of them. At all of these filling frac- 
tions, there exists also one (or several) candidate Abelian 
phase (s); in most cases, it is assumed that these plateaus 
are described by one of the Abc;lian phases. However, 
the pattern-of-zeros construction also yields many sim- 
ple non-Abelian states at these filling fractions. In some 
situations, we expect the non-Abelian states to be good 
candidate states. 

There are at least two dimensionless quantities that 
are important determining factors for which FQH state 
is realized. The first parameter is a = Vlnter/Vintra) where 
Vinter IS the potential for interlayer repulsion and Vlntra 
is the potential for intralayer repulsion. The second pa- 
rameter is 7 = i/Vintra, whcrc t is the interlayer hopping 
amplitude. In the limit a ~ and 7 ~ 0, the system will 
be a FQH state that consists of two independent single- 
layer FQH states in each layer. In the limit 7^1 and 
a ~ 0, a single-layer FQH state may be observed. But if 
we keep 7-^0 and increase a from a ^ 0, then the FQH 
state formed by two independent single- layer FQH states 
in each layer must undergo a phase transition into either a 
compressible phase or a new incompressible state. In the 
latter case, an Abelian hierarchy state (such as a bilayer 
composite fermion state) may form, which would in most 
cases be a state described by a 4 x 4 or more complicated 
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fC-matrix and would have four or more edge modes. The 
other possibihty is that a non-Abehan two-component 
state may form. The pattern-of-zeros construction yields 
non-Abelian two-component wave functions that have ze- 
ros when particles from different layers approach each 
other, indicating that they can accomodate situations in 
which a ~ 1. Additionally, these states generally have 
less than 4 edge modes; if we use the number of edge 
modes as a measure of the complexity of the state, then 
the non-Abelian states are simpler and may therefore be 
realized experimentally. 

At = 2/3, experiments on wide single quantum wells 
have observed a phase transition from a bilaycr to single- 
layer state while experiments on single-layer systems have 
seen a phase transition from a spin-polarized to a spin- 
unpolarized state. In the limit a ~ and 7 ~ 0, the 
system should be in the (3, 3, 0) state. As a is increased 
while 7^0, one possibility is the (1,1,2) state. This 
wave function appears unphysical, because it has higher 
order zeros as particles from different layers approach 
each other than particles from within the same layer. An- 
other wave function, which has the same topological or- 
der as (1, 1, 2), is a spin-singlet composite fcrmion state. '^^ 
There are two other plausible non-Abelian states in this 
situation. One is the following interlayer PfafRan state 
(see (89) and (90) ): 

*2/3|inter= Pf ( ^ ) *(2,2,1) ({^i' ^"0)- (9^) 

The other is the following intralaycr PfafHan state 

(95) 

which has even higher order zeros as particles from dif- 
ferent layers approach each other. ^'2/3|intra has a clus- 
ter structure ^ = 2^ while ^'2/3! inter has a cluster 

structure ^ ~ 2^ ■ At = 2/3 there are also two 

single-layer possibilities that may be realized as 7 is in- 
creased. These arc the particle-hole conjugate of the 
V = XjT) Laughlin state and the Z4 parafermion Read- 
Rezayi state. 



At = 4/5, 4/7, and 4/9, we have the following non- 
Abelian states (see (92) and (93) ): 

*4/5 = ^sc({2:i,'Wi})$(2,2,i)({2;i,f«i}), 
*4/7 = ^sc({2;i,'Wi})*(2,2,|)({2:i,'Wi}), 
*4/9 = $sc({-2i,W'i})*(4,4,i)({-2i.«'i}), (96) 

where = (Ili V'i('2i)V'2(wi)) is a correlation function 
in the .su(3)2/u(l)^ parafermion CFT. These states all 
have 2 1 edge modes. 





Proposed States 


Edge Modes 


Shift 5 




(3,3,0) 


2 


3 




(1,1,2) 


2 


1 


2/3 


2/3|inter (see eqn. (94)) 


2I 


3 




2/3 intra (see eqn. (95)) 


3 


3 




Z4, parafermion 


3 


3 




P-H conjugate oi v = 1/3 


Ih-I-Il 





4/5 


(2/5,2/510) 


4 


4 


s'u(3)2/m(1)^ (see eqn. (96)) 


2| 


3 




(2/3,2/3|l) 


2r + 2l 







(2/7,2/7|0) 


4 


2 


4/7 


su(3)2/m(1)^ (see eqn. (96)) 


2? 


3 


(2/5,2/5|l) 


4 


4 




(2/3, 2/3|2) 


1k-|-3l 







(5,5,3) 


2 


5 


1/4 


(7, 7,1) 

Inter-layer Pfaffian (sec eqn. (98)) 


2 

2| 


7 
7 




Single-layer PfafRan 


1^ 


5 



TABLE I: Proposed explanations for incompressible states 
at experimentally relevant filling fractions, f — 2/3, 4/5, 
4/7, and 1/4, in two-component FQH systems. The bi- 
layer composite fcrmion state {ui, U2\m)'^'^ refers to the state 
Uijizi - ■Wj)"'^i^i{{zi})^L^2{{'Wi}), where is a single layer 
composite fermion state at filling fraction v. For (2/3, 2/3|m), 
we have taken the single layer 2/3 state to be the particle- 
hole conjugate of the Laughlin state, ur + ul indicates that 
there arc tir, right-moving edge modes and ul left-moving 
edge modes. See Appendix C for details of how to calculate 
the number of edge modes and the shift <S. 

The other set of proposed Abclian states arc the bilaycr 
composite fermion states"^"^ (fg, t'olw), which refer to the 
wave function 

*(.o,.o|m) = n(^' - ^jr^'^oi{Zi})^^o{{Wi}). (97) 

Here $^|-, ({zj}) is a single-layer FQH state at filling frac- 
tion 1^0 • These states have 4 edge modes, indicating that 
they may be less stable than the alternative non-Abelian 
possibilities. 

Recently, an incompressible state was found at = 1/4 
and it is unclear what phase this corresponds to and even 
whether it is a single-layer or double-layer phase. Some 
possibilites that have recently been considered^^ are the 
(5,5,3) and (7,7,1) Halpcrin states and the = 1/4 
single-layer PfafHan. The pattern-of-zeros construction 
yields many other alternative possibilities, perhaps the 
most physical (and simplest) of which is the following 
interlayer PfafHan; 

'^{{Zi,Wi})= Pf(-^')$(g_6 2)({^i,t«i}). (98) 

In Table I we summarize some of the filling fractions at 

which incompressible states have been experimentally ob- 
served in two-component FQH systems. For each filling 
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Quasiparticlcs for Intorlayor Pfafiian at v = 2/^ 


i Operator 


iotal Charge 


Scaling Dimension 




1 


3/2 


Ve2 = V'fi V 2"''- 


1 


3/2 


i A / 
3 V 2 ^ — 


2/3 


1/3 


3 V 2 '?*+^^* \/2 ^~ 


1/3 


1/3 


• 1 A3" J, 


1/3 


7/48 




2/3 


31/48 


ae V2 





5/16 







1/2 




2/3 


5/6 




1/3 


5/6 



TABLE II: Quasiparticle operators from the CFT for u = 2/3 
Interlayer Pfaffian states. The scalar boson fields 0+ and <j>- 
arc related to the total and relative density fluctuations of the 
two layers, respectively, a is the spin field in the Ising CFT, 
which has scaling dimension 1/16. 



V 


Charge q^m 


Scaling Dimension h 


2/3 inter 


1/3 


^ + ^ + 


2/.'! iin.ia 


i/(i 


T?i + 4S + TFT 


4/5 


1/5 


- + - + — 

in ^ 40 ^ 94 


4/7 


1/7 


in ^ 5fi ^ s 


4/9 


1/9 


J_ + J. + J_ 

in 72 5R 


1/4 


1/8 





TABLE III: Charge and scaling dimensions of the quasiparti- 
cle operators with minimal nonzero total charge in the non- 
Abclian bilaycr states discussed hero. In the scaling dimen- 
sion, the first term comes from the iion-Abclian part, the 
second term comes from the total density fiuctuations (the 
?7(1) part), the third term comes from the relative density 
fluctuations of the two layers (also the f/(l) part). 



fraction wc list sonic of the proposed wave functions that 
may characterize the topological order of those phases, 
the number of edge modes, and their respective shifts 
on the sphere. We list the qasiparticlcs, their electric 
charges, and their scaUng dimensions for the interlayer 
Pfaffian state at = 2/3 in Table II. In Table III, we 
list the quasiparticles with the minimal electric charge 
and their scaling dimensions h for the non-Abelian FQH 
states discussed in this paper [see eqns. (90), (92), (93), 
(96), (98)]. Those minimally charged quasiparticles may 
dominate interedge tunneling and give rise to the follow- 
ing I-V curve: / a V^'^~^ in the T = limit. 

In summary, we find many simple non-Abelian bilayer 
states that occur at experimentally observed filling frac- 
tions. For certain effective interactions among the elec- 
trons in the unfilled Landau levels, these states may be 
more favorable than their Abelian counterparts. In these 
cases, the non-Abelian states have larger interlayer corre- 
lations and therefore may be energetically more favorable 



f = 2/3 Interlayer Pfaffian 




iVi 

2 4 6 8 10 


2 

6 
8 

10 


(3,3) (7,5) (11,7) (15,9) (19,11) 
(5,7) (9,9) (13,11) (17,13) (21,15) 
(7,11) (11,13) (15,15) (19,17) (23, 19) 
(9,15) (13,17) (17,19) (21,21) (25,23) 
(11,19) (15,21) (19, 23) (23,25) (27,27) 



TABLE IV: Values of (A^^, Ar|) that yield rotationally invari- 
ant states on the sphere for various choices of {Ni,N2) for the 
V = 2/3 interlayer Pfaffian. 



V — 4/5 su(3)2/w(l)^ parafermion state 




iVi 

2 4 6 8 10 


2 

Ar2 ^ 
6 

8 

10 


(2,2) (6,3) (10,4) (14,5) (18,6) 
(3,6) (7,7) (11,8) (15,9) (19,10) 
(4,10) (8,11) (12,12) (6, 13) (20, 14) 
(5,14) (9,15) (13,6) (17,17) (21,18) 
(6,18) (10,19) (14, 20) (18,21) (22,22) 



TABLE V: Values of {N\,,N%) that yield rotationally invari- 
ant states on the sphere for various choices of {Ni,N2) for the 
i> = 4/5 sw(3)2/m(1)^ parafermion state. 



in situations in which the interlayer repulsion is compa- 
rable to the intralayer repulsion. 



A. Conditions on Filling the Sphere 

A useful tool for identifying FQH states in numer- 
ical studies of exact diagonalization on finite systems 
on a sphere is to look at what values of the shift, 
S = y~^Nf. — N^, a ground state with zero total angular 
momentum is found. This then limits the possibilities 
of which topological phase is realized in the system to 
those that have that particular value of the shift. Simi- 
larly, in such numerical studies of multilayer systems, one 
can look for the different sets (A^i, • • • ,Nf; N^, ■ ■ ■ , N^^) 
that yield a ground state with zero total angular mo- 
mentum. Each topological phase will have its own list of 
(iVi, • • • ,Nf;Nl,--- , N^^) that let it fiU the sphere; an- 
alyzing this can be a useful way of determining which 
topological phase is obtained numerically. In Section 
YD, we found conditions that N and iV$ should satisfy 
for the FQH state to fill the sphere. 

^2 o\ 

we find that the condition (46) becomes trivial as long as 
N = 'Y^j kifii, where kj is an integer. This means that 
as long as A^i and N2 are even and N^, satisfy (44), 
then these states can fill the sphere. In this case, we find 



For states that have a cluster structure n = 
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that (44) reduces to the form 

where M is a 2 x 2 matrix and >S is the shift, which 

can be calculated using cqn. (Cll). The states that we 
have been considering are of the form $ = x <&(a,/3,7), 

for which M = {'^ ^ ] . Tables IV and V hsts some 

\(3 a) 

examples. 



XI. SUMMARY 

In this paper, we generalized the pattern-of-zeros char- 
acterization and classification of FQH states to multi- 
component cases. We found that the topological orders 
in a multicomponent FQH state can be characterized by 
the following data: a matrix n that describes the cluster 
structure, a matrix m and a sequence {S^} that describes 
the pattern of zeros. 

Owi pattern-of-zeros characterization gives tis a general 
quantitative view on a large class of Abelian and non- 
Abelian bilayer FQH states, which allow us to determine 
which states are simpler than other states. We find some 
simplest non-Abelian generalizations of the Laughlin- 
Halperin Abelian bilayer states. Those simple non- 
Abelian states may describe some of the bilayer/spin- 
unpolarized FQH states observed in experiments and nu- 
merical calculations. 

This research is partially supported by NSF Grant No. 
DMR-0706078. 



APPENDIX A: OCCUPATION NUMBER 
CHARACTERIZATION 

In the single-component pattern-of-zeros description, 
there is an occupation number characterization that is 
a useful way to understand both the ground states and 
the quasiparticles in FQH states. The generalization to 
multilayer states does not appear to be quite as simple or 
useful, but for the sake of completeness we will analyze 
it below. 

A convenient set of single-particle basis states for par- 
ticles in the lowest Landau level are monomials of the 

form z™, for integer m. Thus, a basis for symmetrized 
wavefunctions of Nj particles of type / is given by: 

N, 

*w}=Enn(4w)'-. (Ai) 

{Pi} I »=1 

where Pj is a permutation of the particles of type 7, 

is an integer, and ni is a vector whose component n\ 
is the number of particles of type / that occupy the l^^ 



orbit. The polynomial that we are interested in can be 
expanded in terms of these basis states as: 

*=EC{na*{na- (A2) 

Now we may ask what kind of boson occupations {n;} 
will be present in the sum (A2) for a polynomial $ with 
a given pattern-of-zeros {5*5}. To answer this question, 
let us set z\ = in <i>({2;/}). Since $ is nonzero when 
z\ = Q due to translation invariance, there must be a 
boson occupation {n;} in the above sum that contains 
at least one boson occupying the (z^)'=° orbital. That 
is, there is a term in the above sum with nj > 0. Now, 
suppose that we bring a second particle of the same type, 
z\ to 0. The minimal power of z\ in <I>(0, • • • ) is Dg^^g^ : 

+ 0((zi)^=-i-i+i). (A3) 

Thus, among those {Hi} that have at least one boson of 
type 1 occupying the (z'^)'^" orbital, there must also be 
an {fti} that contains a second boson of type 1 occupy- 
ing the (z^)'^=i orbital where = Dg^^g^ = S2ei ~ ^Si- 
Next, assume that two bosons occupy the {z^)^ and 
(z^)'2=i orbital, and bring a third particle of type 1 to 
0; the minimal power of Zg is -D2ei,ei: 

P2{zlzl...)r.{zir^^^'^^P,{zlzl...) 

+ 0((4)^-i--i+i). (A4) 

Thus, among those {ni} that have two type 1 bosons oc- 
cupying the I = and ^26*1 orbitals, there is a third boson 
of type 1 occupying the hg^ = £'261, ei = 'S'sei - 'S'2ei o^"- 
bital. Continuing in this way, we see that there must be a 
type 1 boson occupying the orbitals l\g^ = SaSi— S(^a-i)ei 
for a = I,-- - ,A''i. After taking all the type 1 parti- 
cles to 0, we may begin to take the type 2 particles 
to zero, one by one, thus obtaining that there must 
be a type 2 boson occupying the orbitals I'^g^^j^^g^ = 
Sae2+Niei - S(^a-i)e2+Niei- Continuing this argument for 
bosons of every type, we find that there must be a term 
in the sum (A2) with occupation number described by 
the above sequence of /i's. 

However, in the above argument, we chose a particular 
sequence in which to take various particles to zero. We 
first took all of the type 1 particles to zero one-by-one, 
and then all of the type 2 ones, and so on. But we could 
just as well have made the argument with any sequence. 
Suppose that after taking i particles to the origin, there 
is a aj-clustcr at the origin. Thus {di} is a sequence 
that describes the order in which we take particles to 
the origin until all particles are at the origin. For every 
such sequence, we may make the above argument and 
argue that if ai+i = di + ej, then there must be a type / 
boson occupying the orbital ll,^^ = Ssi+Si — Ssi- If we 
enumerate all the different sequences {a^} by an integer 
a, then by considering each a, we see that there must 
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be a term in the sum (A2) with occupation number ft", 
nj'" would be the number of i, along the sequence {oj}, 
for which if = I. Notice that ll must be non-zero. Thus 
we have the following important condition on S3: 



— S3 — Ss-gj > 0. 



(A5) 



The analysis above can be thought of in the following 

way. Consider an A^/-dimensional lattice Z^-f , where Nf 
is the number of layers. At every site a of this lattice 
(a/ > 0) we can associate the nonnegative integer S3. 
On each link (a, a — ej), of the lattice we may also as- 
sociate an integer ll = S3 — Ss-e/- Now consider any 
directed path from the origin to N {Nj is the number of 
particles of type /), in which the sum of the coordinates 
of every point on the path is one larger than the sum of 
the coordinates of the point preceding it, and enumerate 
the set of these paths by a. To each such path we asso- 
ciate an occupation number sequence n", where nj'" is 
the number of links along the path a whoso ll = I. If 
$ has a pattern-of-zeros {S3}, then its basis expansion 
(A2) must contain a term with occupation number n". 
Thus we may rewrite (A2) as 



{n,} 



The two sequences {S3} and {nf} contain the same infor- 
mation and are one-to-one labellings of each other. How- 
ever, {nf} is redundant in the sense that it does not need 
to be specified for every a in order to reconstruct I3. The 
hi that appear in the second sum characterize the sub- 
leading terms that appear when coordinates are brought 
together; thus those {hi} correspond to sequences {S3} 
where S3> S3. 

In the single-layer case, {la} naturally defined an oc- 
cupation number sequence {ni}, which also described 
the FQH state in the thin-cylinder limit. In the multi- 
component generalization, we have {1 3}, which seems to 
admit no simple generalization of the above occupation 
number sequence. Instead, one has such occupation num- 
ber distributions for a large number of sequences which 
we enumerated above by a. We have not analyzed on 
general grounds which particular sequences contribute 
the most weight to the wavefunction in the thin cylin- 
der limit. 



APPENDIX B: su{3)2/u{lf PARAFERMION 
CONFORM AL FIELD THEORY 

Some of the simplest non-Abelian bilayer states are 
closely related to sm(3)2/u(1)^ parafermion CFT. This 
CFT has central charge c — Q/5 and has three simple 
currents, V'a; ^/3) and ipa+f3 = V'aV'/3) all of which square 
to the identity and have scaling dimension 1/2. 

There arc four other primary fields, which arc associ- 
ated with the fundamental representation of su{'S). Their 



Primary fields in su(3)2/it(l)^ 


CFT Operator 


Scaling Dimension 


a 


1/10 




1/10 


afj = tp/ja 


6/10 




1/10 



TABLE VI: Primary fields and their scaling dimensions in the 
sw(3)2/w(l)^ parafermion CFT. 



scaling dimensions are listed in Table VI. The fusion 
rules for these fields all follow from the following fusion 
rule: 



G X a 



1 



■ 17/3. 



(Bl) 



APPENDIX C: CALCULATIONS FOR 
CANDIDATE STATES 

1. Number of Edge Modes 

The total number of edge modes is equal to the cen- 
tral charge of the corresponding CFT for the states that 
are described by the pattern of zeros. For the hierarchy 
states, the number of edge modes is given by the rank of 
the if-matrix. Furthermore, in the latter case, the num- 
ber of right (left) -moving edge modes is given by the 
number positive (negative) eigenvalues of the ii'-matrix. 

The interlayer Pfaffian states are described by a CFT 
that consists of the Ising CFT, with c = 1/2, and two 
scalar boson CFTs, each with c = 1. Thus the number 
of edge modes for the interlayer Pfaffian is 2^. 

The intralayer Pfaffian states have two Ising CFTs in 
addition to the two scalar boson CFTs, so the total num- 
ber of edge modes is 3. 

The central charge of the su(3)2/u(l)^ parafermion 
CFT is c = 6/5; the two-component FQH states based 
on this are described by the su{2>)2/u(Xf theory and two 
scalar bosons, for a total of 2| edge modes. 

For the {m,m,l) states, the if-matrix is K = 

(m I 
I m 

all edge modes move in the same direction; if m < I, 
then there is one right-moving and one left-moving edge 
mode. 

For the states {uq, volm), the K matrix is 



These states have 2 edge modes; if m > /, 



/Al\ K 



K 



'0 

12 



m \ 

Kl, if 
V ifOi K^^J 



^21 ^22 



(CI) 



where is the ff -matrix in the hierarchical basis of the 
state For the v = 2/5 state, AT^ in the hierarchical 
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basis is 



3 -1 
-12, 



(C2) 



For the u = 2/7 state, in the hierarchical basis is 



3 1 ' 

1 -2, 



(C3) 



For the i/ = 2/3 P-H conjugate of the 1/3 Laughhn state, 
in the hierarchical basis is 



1 1 ' 
1 -2, 



(C4) 



Using this, we find that the P-H conjugate of the 1/3 
Laughlin state has edge modes 1^ + 1^, (2/3, 2/3|l) has 
edge modes 2r + 2l, and (2/3, 2/3|2) has 1^ + 3z,. 



2. Shifts on Sphere 

In the hierarchy basis, the formula for the shift is given 
by36 



(C5) 



Using this formula, we find 5 = for the particle-hole 
conjugate of the = 1/3 Laughhn state. 

Now consider the bilayer composite fermion state 
{uo,i'o\m): 



= l[{z, - wjr^^,i{z,})K,{{w,}). (C6) 



Let be the maximum power of Zi in $^,(,({2;,}). It 
satisfies: 



(C7) 



where Sq is the shift of the state The factor 

Hj — Wj)"^ increases the power of zi by mN2. Thus 
the maximum power of Zi in ^(^„,^g\m) is N^: 

Nl = Uo^Ni + mN2 - So . (C8) 



In our cases, Ni = N2, and the number of flux quanta is 
the same in each layer, so 



{lyQ^ + m)Ni - So, 



(C9) 



ly^^ = 1^(7^ + TO is the filling fraction in one layer. Thus 
we see that the shift of ^(^vo^vo\m) is also So- 

For (2/3, 2/3|to), we take $2/3 to be the particle-hole 
conjugate of the 1/3 Laughlin state. $2/3 has shift So = 
0. Thus (2/3,2/3|m) also has shift S = 0. (2/5,2/5|to) 
has shift 4, because $2/5 iias shift 4. 

For the Halperin (m, n, I) states, the iiT-matrix can be 

/to a 

written as K = \ .In this basis, the shift is given 

yi nj 

by 



5 = z.-i^A7;i^^^. 



(CIO) 



IJ 



For the states described by the pattern of zeros, we can 
use the following formula 

^ ^ f i-'^^Y.i VI {mi I - Sfij + Sfii-ei) for bosons 
\ i^~^J2i t^iimii + 1- Sfii + Sfii-ei) for fermions 

(Cll) 



3. Electron and Quasiparticle Operators for 
sw(3)2/w(l)^ states 

The electron operators for the sm(3)2/u(1)^ FQH 
states that we discuss are of the form: 



Ve2 = ^pe'^'>'+-' 



(C12) 



where s = V3/2, 1/2, and ^7/2 for the v = 4/5, 4/7 
and 4/9 states, respectively. The quasiparticle operators 
with minimal total charge are of the form: 



Var 



(C13) 



. The 



and have scaling dimension hqp = + ^ 1 2 
total charge of the quasiparticle is Q = 1/5, 1/7, and 
1/9 for the v = 4/5, 4/7 and 4/9 states, respectively. 
Sqp = 1/3, 1, and 1/7, respectively, for these states. 
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